HMZDupFinder: a robust computational approach for detecting intragenic homozygous duplications from exome sequencing data.

Homozygous duplications contribute to genetic disease by altering gene dosage or disrupting gene regulation and can be more deleterious to organismal biology than heterozygous duplications. Intragenic exonic duplications can result in loss-of-function (LoF) or gain-of-function (GoF) alleles that when homozygosed, i.e. brought to homozygous state at a locus by identity by descent or state, could potentially result in autosomal recessive (AR) rare disease traits. However, the detection and functional interpretation of homozygous duplications from exome sequencing data remains a challenge. We developed a framework algorithm, HMZDupFinder, that is designed to detect exonic homozygous duplications from exome sequencing (ES) data. The HMZDupFinder algorithm can efficiently process large datasets and accurately identifies small intragenic duplications, including those associated with rare disease traits. HMZDupFinder called 965 homozygous duplications with three or less exons from 8,707 ES with a recall rate of 70.9% and a precision of 16.1%. We experimentally confirmed 8/10 rare homozygous duplications. Pathogenicity assessment of these copy number variant alleles allowed clinical genomics contextualization for three homozygous duplications alleles, including two affecting known OMIM disease genes EDAR (MIM# 224900), TNNT1(MIM# 605355), and one variant in a novel candidate disease gene: PAAF1.

Diminished TMEM100 Expression in a Newborn With Acinar Dysplasia and a Novel TBX4 Variant: A Case Report.

Acinar dysplasia (AcDys) of the lung is a rare lethal developmental disorder in neonates characterized by severe respiratory failure and pulmonary arterial hypertension refractory to treatment. Recently, abnormalities of TBX4-FGF10-FGFR2-TMEM100 signaling regulating lung development have been reported in patients with AcDys due to heterozygous single-nucleotide variants or copy-number variant deletions involving TBX4, FGF10, or FGFR2. Here, we describe a female neonate who died at 4 hours of life due to severe respiratory distress related to AcDys diagnosed by postmortem histopathologic evaluation. Genomic analyses revealed a novel deleterious heterozygous missense variant c.728A>C (p.Asn243Thr) in TBX4 that arose de novo on paternal chromosome 17. We also identified 6 candidate hypomorphic rare variants in the TBX4 enhancer in trans to TBX4 coding variant. Gene expression analyses of proband's lung tissue showed a significant reduction of TMEM100 expression with near absence of TMEM100 within the endothelium of arteries and capillaries by immunohistochemistry. These results support the pathogenicity of the detected TBX4 variant and provide further evidence that disrupted signaling between TBX4 and TMEM100 may contribute to severe lung phenotypes in humans, including AcDys.

SNV/indel hypermutator phenotype in biallelic RAD51C variant: Fanconi anemia.

We previously reported a fetus with Fanconi anemia (FA), complementation group O due to compound heterozygous variants involving RAD51C. Interestingly, the trio exome sequencing analysis also detected eight apparent de novo mosaic variants with variant allele fraction (VAF) ranging between 11.5 and 37%. Here, using whole genome sequencing and a 'home-brew' variant filtering pipeline and DeepMosaic module, we investigated the number and signature of de novo heterozygous and mosaic variants and the hypothesis of a rare phenomenon of hypermutation. Eight-hundred-thirty apparent de novo SNVs and 21 de novo indels had VAFs below 37.41% and were considered postzygotic somatic mosaic variants. The VAFs showed a bimodal distribution, with one component having an average VAF of 25% (range: 18.7-37.41%) (n = 446), representing potential postzygotic first mitotic events, and the other component with an average VAF of 12.5% (range 9.55-18.69%) (n = 384), describing potential second mitotic events. No increased rate of CNV formation was observed. The mutational pattern analysis for somatic single base substitution showed SBS40, SBS5, and SBS3 as the top recognized signatures. SBS3 is a known signature associated with homologous recombination-based DNA damage repair error. Our data demonstrate that biallelic RAD51C variants show evidence for defective genomic DNA damage repair and thereby result in a hypermutator phenotype with the accumulation of postzygotic de novo mutations, at least in the prenatal period. This 'genome hypermutator phenomenon' might contribute to the observed hematological manifestations and the predisposition to tumors in patients with FA. We propose that other FA groups should be investigated for genome-wide de novo variants.

SNV/indel hypermutator phenotype in biallelic RAD51C variant - Fanconi anemia.

We previously reported a fetus with Fanconi anemia (FA), complementation group O due to compound heterozygous variants involving RAD51C . Interestingly, the trio exome sequencing analysis also detected eight apparent de novo mosaic variants with variant allele fraction (VAF) ranging between 11.5%-37%. Here, using whole genome sequencing and a 'home-brew' variant filtering pipeline and DeepMosaic module, we investigated the number and signature of de novo heterozygous and mosaic variants and the rare phenomenon of hypermutation. Eight-hundred-thirty apparent dnSNVs and 21 de novo indels had VAFs below 37.41% and were considered postzygotic somatic mosaic variants. The VAFs showed a bimodal distribution, with one component with an average VAF of 25% (range: 18.7-37.41%) (n=446), representing potential postzygotic first mitotic events, and the other component with an average VAF of 12.5% (range: 9.55-18.69%) (n=384), describing potential second mitotic events. No increased rate of CNV formation was observed. The mutational pattern analysis for somatic single base substitution showed SBS40, SBS5, and SBS3 as the top recognized signatures. SBS3 is a known signature associated with homologous recombination-based DNA damage repair error. Our data demonstrate that biallelic RAD51C variants show evidence for defective genomic DNA damage repair and thereby result in a hypermutator phenotype with the accumulation of postzygotic de novo mutations, at least in the prenatal period. This 'genome hypermutator phenomenon' might contribute to the observed hematological manifestations and the predisposition to tumors in patients with FA, and pregnancy loss in general. We propose that other FA groups should be investigated for genome-wide de novo variants.

Coexisting Conditions Modifying Phenotypes of Patients with 22q11.2 Deletion Syndrome.

22q11.2 deletion syndrome (22q11.2DS) is the most common genomic disorder with an extremely broad phenotypic spectrum. The aim of our study was to investigate how often the additional variants in the genome can affect clinical variation among patients with the recurrent deletion. To examine the presence of additional variants affecting the phenotype, we performed microarray in 82 prenatal and 77 postnatal cases and performed exome sequencing in 86 postnatal patients with 22q11.2DS. Within those 159 patients where array was performed, 5 pathogenic and 5 likely pathogenic CNVs were identified outside of the 22q11.2 region. This indicates that in 6.3% cases, additional CNVs most likely contribute to the clinical presentation. Additionally, exome sequencing in 86 patients revealed 3 pathogenic (3.49%) and 5 likely pathogenic (5.81%) SNVs and small CNV. These results show that the extension of diagnostics with genome-wide methods can reveal other clinically relevant changes in patients with 22q11 deletion syndrome.

Cloud-native distributed genomic pileup operations.

MOTIVATION: Pileup analysis is a building block of many bioinformatics pipelines, including variant calling and genotyping. This step tends to become a bottleneck of the entire assay since the straightforward pileup implementations involve processing of all base calls from all alignments sequentially. On the other hand, a distributed version of the algorithm faces the intrinsic challenge of splitting reads-oriented file formats into self-contained partitions to avoid costly data exchange between computational nodes. RESULTS: Here, we present a scalable, distributed and efficient implementation of a pileup algorithm that is suitable for deploying in cloud computing environments. In particular, we implemented: (i) our custom data-partitioning algorithm optimized to work with the alignment reads, (ii) a novel and unique approach to process alignment events from sequencing reads using the MD tags, (iii) the source code micro-optimizations for recurrent operations, and (iv) a modular structure of the algorithm. We have proven that our novel approach consistently and significantly outperforms other state-of-the-art distributed tools in terms of execution time (up to 6.5× faster) and memory usage (up to 2× less), resulting in a substantial cloud cost reduction. SeQuiLa is a cloud-native solution that can be easily deployed using any managed Kubernetes and Hadoop services available in public clouds, like Microsoft Azure Cloud, Google Cloud Platform, or Amazon Web Services. Together with the already implemented distributed range join and coverage calculations, our package provides end-users with a unified SQL interface for convenient analyses of population-scale genomic data in an interactive way. AVAILABILITY AND IMPLEMENTATION: https://biodatageeks.github.io/sequila/.

High-level gonosomal mosaicism for a pathogenic non-coding CNV deletion of the lung-specific FOXF1 enhancer in an unaffected mother of an infant with ACDMPV.

BACKGROUND: Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) results from haploinsufficiency of the mesenchymal transcription factor FOXF1 gene. To date, only one case of an ACDMPV-causative CNV deletion inherited from a very-low level somatic mosaic mother has been reported. METHODS: Clinical, histopathological, and molecular studies, including whole genome sequencing, chromosomal microarray analysis, qPCR, and Sanger sequencing, followed by in vitro fertilization (IVF) with preimplantation genetic testing (PGT) were used to study a family with a deceased neonate with ACDMPV. RESULTS: A pathogenic CNV deletion of the lung-specific FOXF1 enhancer in the proband was found to be inherited from an unaffected mother, 36% mosaic for this deletion in her peripheral blood cells. The qPCR analyses of saliva, buccal cells, urine, nail, and hair samples revealed 19%, 18%, 15%, 19%, and 27% variant allele fraction, respectively, indicating a high recurrence risk. Grandparental studies revealed that the deletion arose on the mother's paternal chromosome 16. PGT studies revealed 44% embryos with the deletion, reflecting high-level germline mosaicism. CONCLUSION: Our data further demonstrate the importance of parental testing in ACDMPV families and reproductive usefulness of IVF with PGT in families with high-level parental gonosomal mosaicism.

Variable expressivity in a four-generation ACDMPV family with a non-coding hypermorphic SNV in trans to the frameshifting FOXF1 variant.

Heterozygous single nucleotide variants (SNVs) or copy-number variant deletions involving FOXF1 or its distant lung-specific enhancer on chromosome 16q24.1 have been identified in 80-90% of patients with Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV), a lethal neonatal lung developmental disorder. We describe a four-generation family with a deceased ACDMPV neonate, her sibling from the electively terminated pregnancy, healthy mother with a history of pulmonary arterial hypertension (PAH), an unaffected aunt, an aunt deceased due to findings consistent with ACDMPV, and a reportedly unaffected grandmother, all with the frameshifting variant c.881_902dup (p.Gly302Profs*46) in FOXF1, and a deceased great-grandmother with a history of PAH. Genome sequencing analyses in the proband's unaffected mother revealed a non-coding putative regulatory SNV rs560517434-A within the lung-specific distant FOXF1 enhancer in trans to the FOXF1 frameshift mutation. Functional testing of this variant using an in vitro luciferase reporter assay showed that it increased FOXF1 promoter activity 10-fold. Our studies further demonstrate that non-coding SNVs in the FOXF1 enhancer region can rescue the lethal ACDMPV phenotype and support the compound inheritance gene dosage model.

Loss-of-function variant in chymotrypsin like elastase 3B (CELA3B) is associated with non-alcoholic chronic pancreatitis.

BACKGROUND: Genetic alterations in digestive enzymes have been associated with chronic pancreatitis (CP). Recently, chymotrypsin like elastase 3B (CELA3B) emerged as a novel risk gene. Thus, we evaluated CELA3B in two European cohorts with CP. METHODS: We analyzed all 8 CELA3B exons in 550 German non-alcoholic CP (NACP) patients and in 241 German controls by targeted DNA sequencing. In addition, we analyzed exons 6 and 7 by Sanger sequencing and the c.129+1G>A variant by melting curve analysis in 1078 further German controls. As replication cohort, we investigated up to 243 non-German European NACP patients and up to 1665 controls originating from Poland, Hungary, and Sweden. We assessed the cellular secretion and the elastase activity of recombinant CELA3B variants. RESULTS: In the German discovery cohort, we detected a splice-site variant in intron 2, c.129+1G>A, in 9/550 (1.64%) CP patients and in 5/1319 (0.38%) controls (P=0.007, OR=4.4, 95% CI=1.5-13.0). In the European replication cohort, this variant was also enriched in patients (9/178 [5.06%]) versus controls (13/1247 [1.04%]) (P=0.001, OR=5.1, 95% CI=2.1-12.0). We did not find the two previously reported codon 90 variants, p.R90C and p.R90L. CONCLUSIONS: Our data indicate that CELA3B is a susceptibility gene for CP. In contrast to previous reports suggesting that increased CELA3B activity is associated with CP risk, the splice-site variant identified here is predicted to cause diminished CELA3B expression. How reduced CELA3B function predisposes to pancreatitis remains to be elucidated.

Postzygotic mutations and where to find them - Recent advances and future implications in the field of non-neoplastic somatic mosaicism.

The technological progress of massively parallel sequencing (MPS) has triggered a remarkable development in the research on postzygotic mutations. Although the overwhelming majority of studies in the field focus on oncogenesis, non-neoplastic diseases are attracting more and more attention. The aim of this review was to summarize some of the most recent findings in the field of somatic mosaicism in diseases other than neoplastic events. We discuss the abundance and role of postzygotic mutations, with a special emphasis on disorders which occur only in a mosaic form (obligatory mosaic diseases; OMDs). Based on the list of OMDs compiled from the published literature and three databases (OMIM, Orphanet and MosaicBase), we demonstrate the prevalence of cancer-related genes across OMDs and suggest other sources to further explore OMDs and OMD-related genes. Additionally, we comment on some practical aspects related to mosaic diseases, such as approaches to tissue sampling, the MPS coverage required to detect variants at a very low frequency, as well as on bioinformatic and molecular tools dedicated to detect somatic mutations in MPS data.

The Thousand Polish Genomes-A Database of Polish Variant Allele Frequencies.

Although Slavic populations account for over 4.5% of world inhabitants, no centralised, open-source reference database of genetic variation of any Slavic population exists to date. Such data are crucial for clinical genetics, biomedical research, as well as archeological and historical studies. The Polish population, which is homogenous and sedentary in its nature but influenced by many migrations of the past, is unique and could serve as a genetic reference for the Slavic nations. In this study, we analysed whole genomes of 1222 Poles to identify and genotype a wide spectrum of genomic variation, such as small and structural variants, runs of homozygosity, mitochondrial haplogroups, and de novo variants. Common variant analyses showed that the Polish cohort is highly homogenous and shares ancestry with other European populations. In rare variant analyses, we identified 32 autosomal-recessive genes with significantly different frequencies of pathogenic alleles in the Polish population as compared to the non-Finish Europeans, including C2, TGM5, NUP93, C19orf12, and PROP1. The allele frequencies for small and structural variants, calculated for 1076 unrelated individuals, are released publicly as The Thousand Polish Genomes database, and will contribute to the worldwide genomic resources available to researchers and clinicians.

Do paternal deletions involving the FOXF1 locus on chromosome 16q24.1 manifest with more severe non-lung anomalies?

Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) is a rare lethal lung developmental disorder in neonates due to heterozygous loss-of-function of the mesenchymal transcription factor gene, FOXF1. Interestingly, unlike ACDMPV-causing point mutations in FOXF1 that can be inherited from the mother or father, causative copy-number variant (CNV) deletions arise de novo and almost exclusively on chromosome 16 inherited from the mother (n = 50 vs. n = 3). Here, we describe a fourth case of a de novo paternal CNV deletion encompassing FOXF1, its neighboring long non-coding RNA gene FENDRR, and their distant lung-specific enhancer, identified in a 21-week-old fetus with tetralogy of Fallot, gastrointestinal and genitourinary abnormalities, a single umbilical artery, and patchy histopathological findings of ACDMPV in autopsy lung. We also review the ACDMPV-causative CNV deletions detected prenatally and propose that the majority of paternal deletions manifest with more severe additional non-lung abnormalities.

Variants in the pancreatic CUB and zona pellucida-like domains 1 (CUZD1) gene in early-onset chronic pancreatitis - A possible new susceptibility gene.

OBJECTIVE: Non-alcoholic chronic pancreatitis (NACP) frequently develops in the setting of genetic susceptibility associated with alterations in genes that are highly expressed in the pancreas. However, the genetic basis of NACP remains unresolved in a significant number of patients warranting a search for further risk genes. DESIGN: We analyzed CUZD1, which encodes the CUB and zona pellucida-like domains 1 protein that is found in high levels in pancreatic acinar cells. We sequenced the coding region in 1163 European patients and 2018 European controls. In addition, we analyzed 297 patients and 1070 controls from Japan. We analyzed secretion of wild-type and mutant CUZD1 from transfected cells using Western blotting. RESULTS: In the European cohort, we detected 30 non-synonymous variants. Using different prediction tools (SIFT, CADD, PROVEAN, PredictSNP) or the combination of these tools, we found accumulation of predicted deleterious variants in patients (p-value range 0.002-0.013; OR range 3.1-5.2). No association was found in the Japanese cohort, in which 13 non-synonymous variants were detected. Functional studies revealed >50% reduced secretion of 7 variants, however, these variants were not significantly enriched in European CP patients. CONCLUSION: Our data indicate that CUZD1 might be a novel susceptibility gene for NACP. How these variants predispose to pancreatitis remains to be elucidated.

Exacerbation of mild lung disorders to lethal pulmonary hypoplasia by a noncoding hypomorphic SNV in a lung-specific enhancer in trans to the frameshifting TBX4 variant.

Variants involving TBX4 are associated with a wide variety of disorders, including pulmonary arterial hypertension, ischiocoxopodopatellar syndrome (ICPPS)/small patella syndrome (SPS), lethal lung developmental disorders (LLDDs) in neonates, heart defects, and prenatally lethal posterior amelia with pelvic and pulmonary hypoplasia syndrome. The objective of our study was to elucidate the wide variable phenotypic expressivity and incomplete penetrance in a three-generation family with a truncating variant in TBX4. In addition to exome and genome sequencing analyses, a candidate noncoding regulatory single nucleotide variant (SNV) within the lung-specific TBX4 enhancer was functionally tested using an in vitro luciferase reporter assay. A heterozygous frameshift variant c.1112dup (p.Pro372Serfs*14) in TBX4 was identified in patients with mild interstitial lung disease (1), bronchiolitis obliterans (1), recurrent pneumothorax (1), ICPPS/SPS (1), LLDD (2), and in unaffected individuals (4). In two deceased neonates with LLDD, we identified a noncoding SNV rs62069651-C located in trans to the mutated TBX4 allele that reduced the TBX4 promoter activity by 63% in the reporter assay. Our findings provide a functional evidence for the recently reported model of complex compound inheritance in which both TBX4 coding and in trans noncoding hypomorphic variants in the lung-specific enhancer of TBX4 contribute to LLDD.

Detection of low-level parental somatic mosaicism for clinically relevant SNVs and indels identified in a large exome sequencing dataset.

BACKGROUND: Due to the limitations of the current routine diagnostic methods, low-level somatic mosaicism with variant allele fraction (VAF) < 10% is often undetected in clinical settings. To date, only a few studies have attempted to analyze tissue distribution of low-level parental mosaicism in a large clinical exome sequencing (ES) cohort. METHODS: Using a customized bioinformatics pipeline, we analyzed apparent de novo single-nucleotide variants or indels identified in the affected probands in ES trio data at Baylor Genetics clinical laboratories. Clinically relevant variants with VAFs between 30 and 70% in probands and lower than 10% in one parent were studied. DNA samples extracted from saliva, buccal cells, redrawn peripheral blood, urine, hair follicles, and nail, representing all three germ layers, were tested using PCR amplicon next-generation sequencing (amplicon NGS) and droplet digital PCR (ddPCR). RESULTS: In a cohort of 592 clinical ES trios, we found 61 trios, each with one parent suspected of low-level mosaicism. In 21 parents, the variants were validated using amplicon NGS and seven of them by ddPCR in peripheral blood DNA samples. The parental VAFs in blood samples varied between 0.08 and 9%. The distribution of VAFs in additional tissues ranged from 0.03% in hair follicles to 9% in re-drawn peripheral blood. CONCLUSIONS: Our study illustrates the importance of analyzing ES data using sensitive computational and molecular methods for low-level parental somatic mosaicism for clinically relevant variants previously diagnosed in routine clinical diagnostics as apparent de novo.

Exome Sequencing Reveals Novel Variants and Expands the Genetic Landscape for Congenital Microcephaly.

Congenital microcephaly causes smaller than average head circumference relative to age, sex and ethnicity and is most usually associated with a variety of neurodevelopmental disorders. The underlying etiology is highly heterogeneous and can be either environmental or genetic. Disruption of any one of multiple biological processes, such as those underlying neurogenesis, cell cycle and division, DNA repair or transcription regulation, can result in microcephaly. This etiological heterogeneity manifests in a clinical variability and presents a major diagnostic and therapeutic challenge, leaving an unacceptably large proportion of over half of microcephaly patients without molecular diagnosis. To elucidate the clinical and genetic landscapes of congenital microcephaly, we sequenced the exomes of 191 clinically diagnosed patients with microcephaly as one of the features. We established a molecular basis for microcephaly in 71 patients (37%), and detected novel variants in five high confidence candidate genes previously unassociated with this condition. We report a large number of patients with mutations in tubulin-related genes in our cohort as well as higher incidence of pathogenic mutations in MCPH genes. Our study expands the phenotypic and genetic landscape of microcephaly, facilitating differential clinical diagnoses for disorders associated with most commonly disrupted genes in our cohort.

Loss of function TRPV6 variants are associated with chronic pancreatitis in nonalcoholic early-onset Polish and German patients.

PURPOSE: Loss of function variants of the transient receptor potential cation channel, subfamily V, member 6 (TRPV6) have been recently associated with chronic pancreatitis (CP) in Japanese, German and French patients. Here, we investigated the association of TRPV6 variants with CP in independent European cohorts of early-onset CP patients from Poland and Germany. PATIENTS AND METHODS: We enrolled 152 pediatric CP patients (median age 8.6 yrs) with no history of alcohol/smoking abuse and 472 controls from Poland as well as 157 nonalcoholic young CP patients (median age 20 yrs) and 750 controls from Germany. Coding regions of TRPV6 were screened by Sanger and next generation sequencing. Selected, potentially pathogenic TRPV6 variants were expressed in HEK293T cells and TRPV6 activity was analyzed using ratiometric Ca2+ measurements. RESULTS: Overall, we identified 10 novel (3 nonsense and 7 missenses) TRPV6 variants in CP patients. TRPV6 p.V239SfsX53 nonsense variant and the variants showing significant decrease in intracellular Ca2+ concentration in HEK293T cells (p.R174X, p.L576R, p.R342Q), were significantly overrepresented in Polish patients as compared to controls (6/152, 3.9% vs. 0/358, 0%; P = 0,0007). Nonsense TRPV6 variants predicted as loss of function (p.V239SfsX53 and p.R624X) were also significantly overrepresented in German patients (3/157; 2.0% vs 0/750; 0%, P = 0.005). CONCLUSIONS: We showed that TRPV6 loss of function variants are associated with elevated CP risk in early-onset Polish and German patients confirming that TRPV6 is a novel CP susceptibility gene.

Wide Fontanels, Delayed Speech Development and Hoarse Voice as Useful Signs in the Diagnosis of KBG Syndrome: A Clinical Description of 23 Cases with Pathogenic Variants Involving the ANKRD11 Gene or Submicroscopic Chromosomal Rearrangements of 16q24.3.

KBG syndrome is a neurodevelopmental autosomal dominant disorder characterized by short stature, macrodontia, developmental delay, behavioral problems, speech delay and delayed closing of fontanels. Most patients with KBG syndrome are found to have a mutation in the ANKRD11 gene or a chromosomal rearrangement involving this gene. We hereby present clinical evaluations of 23 patients aged 4 months to 26 years manifesting clinical features of KBG syndrome. Mutation analysis in the patients was performed using panel or exome sequencing and array CGH. Besides possessing dysmorphic features typical of the KBG syndrome, nearly all patients had psychomotor hyperactivity (86%), 81% had delayed speech, 61% had poor weight gain, 56% had delayed closure of fontanel and 56% had a hoarse voice. Macrodontia and a height range of -1 SDs to -2 SDs were noted in about half of the patients; only two patients presented with short stature below -3 SDs. The fact that wide, delayed closing fontanels were observed in more than half of our patients with KBG syndrome confirms the role of the ANKRD11 gene in skull formation and suture fusion. This clinical feature could be key to the diagnosis of KBG syndrome, especially in young children. Hoarse voice is a previously undescribed phenotype of KBG syndrome and could further reinforce clinical diagnosis.